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We determine diffractive parton distributions of the proton from DGLAP based fits 
to HERA data including the twist-4 contribution from longitudinal polarized virtual 
photons, which is known to be important in the region of large (3. The biggest im- 
pact of this contribution is on the diffractive gluon distribution and on the diffractive 
longitudinal structure function F£ to be determined from HERA data. 



1 Introduction 

Diffractive deep inelastic scattering (DDIS) at HERA provides a very interesting example 
of processes with a clear interplay between hard and soft aspects of the QCD interactions. 
In these interactions, the scattered proton remains intact and a diffractive system is formed 
which is separated in rapidity from the scattered proton. The diffractive events in DIS are 
not rare, on the contrary, they form around 10% of all DIS events [11[31I31[S]. 

After the integration over the proton azimuthal angle, the diffractive cross section is 
given in terms of two structure functions, and Fj^ , which depend on four variables: 
Bjorken-x, photon virtuality and 

+ ' (1) 

Here M is mass of the diffractive system, W is invariant energy of the gamma-proton 
system and p, p' are incident and scattered proton momenta, respectively. In our analysis, 
the diffractive structure functions arc decomposed into the leading twist-2 and higher twist 
contributions 

F2^L - <L + F^^'""^^ + .... (2) 

In the Bjorkcn limit, the leading part depends logarithmically on while the twist-4 part 
is suppressed by an additional power of l/Q^. However, this contribution dominates over 
the twist-2 one for small diffractive masses, M"^ <^ , playing especially important role 
in DDIS. Physically, the twist-4 contribution is given by diffractively produced qq pairs 
from longitudinally polarized virtual photons. The effect of this contribution is particularly 
important for the diffractive longitudinal structure function F^ which is supposed to be 
determined from the HERA data. In [5] we have studied in detail the impact of the higher 
twist contribution on the determination of diffractive parton distributions and the structure 
function F^ . For an alternative non-standard analysis of diffractive DIS, see [7]. 
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2 Twist-2 contribution 



The twist-2 contribution is given in terms of the difFractive parton distributions (DPD) 
through standard coUinear factorization formulae [H [9l [101 E] • For example, in the leading 
logarithmic approximation 

Nf 

pDit^^) =J2e}p{q{,ix,Q^Xp,t) + q{,ix,Q',Xjp,t)} (3) 
/=i 

and Ff = 0. Fr om symmetry of the vacuum quantum number exchange which leads 
to diffraction, we have q(j = qjj = J^/{2Nf) where S is a singlet quark distribution. In 
addition, we also consider a diffractivc gluon distribution g^ which enters the formulae for 
the structure functions in the next-to-leading logarithmic approximation. 

The evolution of the DPD with is governed by the DGLAP equations, and we assume 
Regge factorization of the DPD 

S = fjp{xp,t)^p{p,Q^) (4) 
9D = fp{xp,t) gp{P,Q'^) (5) 

where /3 = x/xp plays the role of the Bjorken variable for DIS diffraction. The motivation 
for such a factorization is a model of difFractive interactions with the pomeron exchange 
|12j . In this model fp is the pomeron flux 

Jip(xp,t) ^ N -^-^Xjp (6) 

where ap{t) = ap(J)) + a'pt is the pomeron Regge trajectory and Fp{t) = Fp{0) e~^o'*l is 
the elastic formfactor which describes the pomeron-proton coupling. The difFractive slope is 
taken form HERA data, Bd = 5.5 GeV~^ and F^{0) = 54.4 GeV"^ [10]. For the pomeron 
trajectory, a'p = 0.25 GeV~^ but the intercept ap{0) is fitted to the analyzed data. 

The so-called pomeron parton distributions, "Sp and gp, are fitted to HERA data as- 
suming the following form of these distributions at the initial scale Qq = 1.5 GeV: 

EpiP) - Ag p""--' {I ~ (3f- (7) 
gp{P) ^ Ag (3^^-^ {I - (3f^ . (8) 

In our fits we used the next-to- leading order DGLAP equations with the value of Aqcd = 
407 MeV (for Nf = 3 flavours) in the strong coupling constant. We also included the cc 
production from the photon-gluon fusion process. In the considered approximation, the 
diffractive structure function pj^^*"^'^^ ^ and it is given as a convolution of the DPD with 
appropriate coefficient functions [5]. 

We performed fits to HI and ZEUS data separately. For each data set we considered two 
scenarios: with and without the higher twist contribution. 

3 Higher twist contribution 

Analyzing the diffractive DIS by considering successive Fock components of the diffractive 
states, qq, qqg from longitudinal and transversely polarized photons, it was found that 
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for /3 « 1 (i.e. small difFractive masses, Af^ ^ Q^) the higher twist Lqq contribution 
from longitudinal photons dominates over the leading twist-2 terms from transverse photons 
[131 [Ml [15] . This indicates that contrary to the fully inclusive case, the higher twist Lqq 
component cannot be neglected in the region of large (3 in the DGLAP based analysis of the 
previous section. Thus we add the following formula to the standard twist-2 relations [6] 

0^(1-/3) 
475 



4/3 P 



i-/3g2 



with 



where i^o and Jo are Bessel functions. Strictly speaking, formula ^ contains all powers 
of but the twist-4 part, proportional to dominates. The dipole cross section 

a{x]p, r) in cq. (jlOp describes the diffractive interaction of the qq pair, treated as a colorless 
dipole, with the proton. We use the phenomenological successful parameterization of a from 
[TBI for this interaction. 



4 Fit results 



The DPD obtained from fits to HI data in the two scenarios: with and without higher twist 
component, arc shown in Fig. [TJ Wc plot the distributions /3Ejp(/3, Q^) and (3g]p{(3,Q'^) 
for several values of Q^. We see that the quark distributions are practically the same 
while the gluon distribution from the fit with higher twist is strongly peaked near (3 — 1. 



This somewhat surprising result can be 
understood by analyzing the logarithmic 
slope of for fixed /3. From the DGLAP 
equations, we schematically have 

QpD r 
-g^^ ~ Pqq®^IP + PqG^Gp J Pqq 

(11) 

where the negative term sums virtual cor- 
rections. For large f3, the measured slope 
is negative which means that the negative 
term in cq. (|lip must dominate over the 
positive ones. The addition of the higher 
twist contribution to , proportional to 
powers of 1/Q^, contributes negative value 
to the slope. This has to be compensated 
by a larger gluon distribution near /3 = 1 in 
the second term on the r.h.s. of eq. pT|) in 
order to describe the same data. 
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Figure 1: Diffractive parton distributions from 
fits to HI data. 
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In Fig. [2] we show the longitudinal 
diffractive structure function Fj^ from the 
fits with (dashed lines) and without (solid 
lines) higher twist. The upper lines in each 
set of curves correspond to the fits to the 
HI data while the lower ones to the ZEUS 
data. We sec that the contribution from the 
higher twist Lqq component is very impor- 
tant for the prediction of the form of Ff' 
for large /3 (dashed lines). In contrast to 
the pure DGLAP analysis (solid lines), Fj^ 
is not negligible when /3 — > 1. 

5 Summary 

We performed the analysis of the diffractive Figure 2: Predictions for from fits to 
data from HERA. In addition to the stan- HERA data, 
dard twist-2 formulae, we also considered 
higher twist contribution produced by lon- 
gitudinally polarized virtual photons. This 

contribution leads to the diffractive gluon distribution which is stronger peaked near (3=1 
than in the pure leading twist analysis. Moreover, the higher twist significantly enhances 
Ff in the region of /3 > 0.7. This is the most important result of the performed analysis 
which hopefully will be tested against the HERA data. 
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